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(N ■ Abstract 

^ ■ Let {E a {N)} a gz N be a family of |Sjv| = 2 N centered unit Gaussian 

o ' 

^ ■ c N (a,r) := Av (E a (N)E T (N)), and H N (a) = —\^NE a (N) the corre- 

sponding random Hamiltonian. Then the quenched thermodynamical 
Q^' limit exists if, for every decomposition N = Ni + N2, and all pairs 

(<j,t) GSjvX Sat: 



random variables defined by the covariance matrix Cat of elements 



iVi N 2 
c N (v,t) < — CN^TT^a),^)) + — C7v 2 (7r 2 (o-),7r 2 (r)) 



^ ' where ftki?)-, k = 1, 2 are the projections of a G Sjv into £/v fc - The con- 

dition is explicitly verified for the Sherrington-Kirckpatrick, the even 
p-spin, the Derrida REM and the Derrida-Gardner GREM models. 



1 Introduction, Definitions and Results 

It has recently been proved by Guerra and Toninelli ||GuTo|| that for the 
Sherrington-Kirckpatrick (hereafter SK) model (as well as for the even-p- 
spin models) the thermodynamical limit exists for the quenched free energy 
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and almost everywhere for its random realizations. In this paper we single 
out general sufficient conditions that imply the existence of the quenched 
thermodynamical limit for any correlated Gaussian random energy model. 
Our analysis thus includes as special cases not only the even p spin models 



(in particular the SK one, p — 2) but also the Derrida REM model fDel|j , [ De2 
and the Derrida-Gardner GREM fDeG"a 



The paper is organized as follows: in this section we introduce the defini- 
tions and state the results. In section 3, after introducing and elucidating the 
operation of lifting for a family of Gaussian random variables, we describe 
the proof of our theorem. In section 4 we show how our analysis applies to 
the specific examples listed above. 

To define the set up we consider a disordered model having 2 N energy 
levels where N is the size of the system. We label the energy levels by the 
index a = {a±, o~2, • ■ ■ , 0"tv} where each G{ takes the values ±1 for i = 1, . . . , N. 
We denote S^r the set of all a. Then |Sjv| = 2 . Clearly coincides with 
the space of all possible 2 N Ising configurations of length N. 

Definition 1 Denote {E a (N)} a£ -z N a family of 2 N centered unit Gaussian 
random variables: 

Av(£ CT (A0) = 0, (1) 
and covariance matrix Cn with elements defined by 

c N (a,a):=Av(El(N))=l, (2) 

c N (a,r) :=Av(E a (N)E T (N)) . (3) 

Here Av (— ) denotes expectation with respect to the probability measure 

dP(E 1 ,...,E 2N ) = 1 e -±(E,c-iE) dE dE (4) 

1 ^ ^{2ny N det{C) K 1 



Definition 2 

1. For each N the Hamiltonian is given by 

H N (a) = -y/NE^N) . (5) 

2. The partition function of the system is: 

Z N (f3, E) = J2 ^ HN{a) = E e*^W (6) 

3. The quenched free energy /jv(/3) of the system is defined as: 

- Pf N (f3) := a N (/3) := 1 Av (lnZ N ((3, E)) . (7) 

Remark 1 From now on we write E a (N) = E a , dropping the iV-dependence. 
Remark moreover that Definition 1 includes Gaussian families of the form 

E a (N) = J + E ^ i(7i + E ^ i ^ ai °j + E J'hhkViV'jVk + 
i i,j i,j,k 

+ • • • + ^ ] Ji 1 ,i2,...,i N °~ii ~i2--- a iN (^) 

»1,*2, — >»JV 

in which every J is an indipendent Gaussian variable. 
Examples. 

1. The SK model. Consider first the model defined by 

1 N 

E ° := n E Ji J°" i(T i ( 9 ) 

where the Jjj are N 2 i.i.d. unit Gaussian random variables. A short 
computation yields 

Av(E a E T ) = [q N (a,r)] 2 



where, as usual 

1 N 

q N (v,T) := —^2(T k r k (10) 

k=l 

is the overlap between the a and r spin configurations. The standard 
SK model is instead defined by 

1 N 

i<j=i 

However the quenched free energy densities (|7]) of the two models co- 
incide up to a rescaling of the temperature, i.e.: 



a N 



(V2/3) = a N {(3) , (12) 



In fact, Ji^OiOj are centered, unit and i.i.d. Gaussian random variables 
V(i,j), and Jija^aj = Jj t i<jj(Ti. Hence Jija^aj + JjjGjUi = V2Ji t j(Ti(Tj 
(here = denotes equality in distribution of two random variables). 
Therefore, taking into account also the N diagonal terms: 

y/NE a = VnV2E* k + J , (13) 



where J is a centered unit Gaussian variable. By ((^,0) formula fll3l) 
immediately yields the relation (|12|). 



2. The p-spin models. Here we consider the model: 

rr - 

E ° := \lW ^ Jhr-A^h---^ ( 14 ) 

ii,...,» p =l 

where the are once more i.i.d. unit Gaussian random variables. 

As before, a short computation yields 

Av(E CT E r ) = [q N (a,r)] p (15) 



3. The Derrida REM. Here the model is specified by Definition |I]with 

Av(E a E T )=6(a,r) (16) 



4. The Derrida- Gardner GREM. Its inclusion into the above framework 
is described in detail in Section 3.3. 

Definition 3 For each a G X^y letni and7t 2 be the two canonical projections 
over the two subsets X n ± and X n 2 , generated by a partition V of the coordi- 
nates (<7i, . . . , ctat) into a subset of N\ coordinates and into a complementary 
set of N 2 coordinates: Ni + N 2 = N, T, N = T, Nl x T, N2 , tti <S> tt 2 = ls N - 

(Example: Af = 4; a 6 E4 with coordinates denoted {o~\, a 2 , 03, 4 }. Consider 
for Ni = N 2 = 2 the partition Va = (a\, a 2 )U(a 3 , cr 4 ). Then X/y = X^ x Xjv 2 
and the two projections 7Tfc : Xjv — > Xjv fe , k — 1, 2 act in the following way: 
71-1(0-1,02, 0-3,0-4) = (0-1,0-2) and 7^(01,02,03,04) = (0-3,04))- 
Our main result is the following: 

THEOREM 1 Let the covariance matrices Cn fulfill the condition: 

Ni N 2 
c N (a,r) - — c 7Vl (7Ti(0),7ri(r)) - — c Na {ir 2 {<r), 7T 2 (r)) < 0, (17) 

for every N > N, every (a, r) G X^r x Xjv o,nd every decomposition Ni + N 2 = 
N. Then the thermodynamical limit exists, in the sense that 

lim ±Av{logZ N (J3)) = sup ±Av(log Z N (J3)) . (18) 



Remark 2 The result ([lq) can be extended to the almost-everywhere con- 
vergence of free energy density, internal energy and ground state energy with 
elementary probability methods (see | GuTc|| ) 



Remark 3 The conditions ( |I7D are not necessary. The proof itself will show 
that we only need the sign of the quantity in the left hand side of flT^) in 
average, not pointwise. Moreover the condition (|]) can be replaced by a more 
general small deviation vanishing for large N and (0) by a uniform (in N) 
bound over the diagonal terms. We plan to return over such a general case 
elsewhere. 

Remark 4 It is still an open interesting question whether the class of models 
we control the thermodynamical limit of do have, in that limit, the properties 



axiomatically introduced by Ruelle in |[Rul|| to define directly the infinite 



particle systems. To this purpose see ||BS|| , [ [BoKul|| , ||BoKu2|| and [[BoKu3|| 



2 Proof 

Within this section it is useful to consider 2 identical copies of the same 
system: the system 1 is assigned the Hamiltonian H(o~) and the system 2 the 
Hamiltonian H(r). 



Definition 4 The quenched measure over the two copies (—) is defined by 
<-> = Av[Z{p, E)]- 2 ~ S H W+ H W . (19) 

(o-,r)eSjvxS JV 

The definition may of course be generalized to r copies. 

We want now to embed a Gaussian system {E a }^ K into a larger one 
{E T }z L for some K < L. In particular we want to embed two of them of 
size N\ and N 2 into one of size N = Aq + N 2 . Our embedding procedure is 
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defined in terms of the two canonical projections TTj, j = 1,2 from Ejv to Sjy- 
given m Definition (|). 

Definition 5 Given the family {E^}^ of size N\ we lift it to one of size 
N: {E£'}x N defining 

E$ = E M „) • (20) 

Moreover starting from {E ij }y, N2 we define in the same way {E a }z N by 

Ef ] = E Ma) . (21) 



Having defined each family {E a }j] N , {E^}% N and {E& }t; N2 we specify thei 
joint distribution requiring mutual independence. 

Remark 5 : The embedded Gaussian systems {E&}e n and {E^}y: n are 
degenerate: In fact for all a and r such that 7ri(er) = 7Pl(t) 

= E<p . (22) 

Summarizing we define the joint measure of {E a }-£ N , {Ea }s N and {Ea}s N 
dP = dPdP\dP2 defined by the three covariances C^v, Cjvi and Cn 2 - 

Proof of THEOREM 1. 

We proceed in three lemmas. 
Lemma 0. Interpolation 

Given a pair (71*1,^2) as before, following [puTof , we pick three independent 



eir 



Gaussian systems E^ a y j = 0, 1, 2 and introduce the quantity (7r (o~) = a) 

2 

H {NjNltN2) (a,t) := -J2VWjEI% } (23) 

j=0 
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where t = t and t\ = t 2 = (1 — t), and the correspondent partition sum 

Z N (t,/3) := e'^^^. (24) 

It is now easy to see that: 

Z N (1,P) = Z N (P), (25) 

and 

Z N (0,p) = > e ,ri(<T) 7r2(CT) 



tSSjv 2 creSjv; 7T2(ir)=r 

= Z Nl (P)-Z N M (26) 

Lemma 1. Boundedness 

The Jensen inequality 

Av(logZ) < log(Av(Z)) (27) 

implies 

^Av (log Z N (/3)) < log(2) + ^- (28) 

because by (^j) Av (Z) = 2e^ 2 ^ 2 after performing the Gaussian integration. 
Lemma 2. Monotonicity 

Taking the t derivative of the logarithm of (p4]) we get: (here we abbreviate 
Hn,N!,n 2 — H) 

- ^ £ (t.^e— ) , (29) 
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where e = 1 and e\ = e 2 = — 1. 

We now use the integration by parts formula for correlated Gaussian 
variables with covariance Cjj, which states 

df 



Avfa.f) = Av [J2 c j,K 



This yields 



7 CT eS]v fc=o v ft \ iV w 

2 fW I d e~ l3H 

= EE efc \/~r Av E c N k ( n k(v),T k ) ■ 



(tSSjv k=0 \ r fe eSjv fe 

Given now G Ejv fe fixed, we calculate 



dE M Z N (t) 



(30) 



(31) 



d e-W 
dE%> Z N (t) 



a,t) VNkTk ■ Z N (t) - e-W-.*) ■ ^ 



/3- 



9Z N 
8E&> 



z 2 N (t) 



p 



5eS]v,7r fe (^)=T fc 



z N {t) 



The term with A; = in formula ([31]) is easy to calculate and we get: 



N ^ Av ( E E c ^ r ) 



55- 



J N 



-P(H(£,t)+H(a,t)) 



cat (a, a) 

iV/3(l - Ctf(<7,T))t , 



-PH(<r,t) 



^2 c N (a,r)e 

(a",r)eEjvxSjv 



■/3(H(r,t)+H(a,t)) 

(32) 



where < — >t is the quenched measure with respect to the Hamiltonian (p3[) . 
In the same way for the term k = 1 (and similarly for k = 2) we obtain: 
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5L 



-/3H(a,t) 



iVi(l - c iVl (7r 1 (a),7r 1 (r))) t . 



/3{H{i,t)+H(a,t)) 



Summing up the three contributions we obtain: 



~Av( W) ) 



/3 < c N (a,r) - —c Ni (tti (0-),7Ti(t)) -— c N2 (ir 2 (cr), 7r 2 (r)) > t , (34) 



iV 



and, by the hypothesis (|T^): 



Av (log Z N (t)) >0. 



(35) 



Formula (^) together with the boundary conditions 
every JVi + N 2 = N 



and (|26|) gives for 
(36) 



This entails Theorem [1] as explained for instance in [ Hu2 



Remark 6 Lemma 3 is indeed a particular case of a theorem by J-P. Kahane 



(see also |[LT|| , Theorem 3.11, p. 74). The Gaussian process X of |KJ can 
in fact be identified with our Gaussian process y/~N E, and the process Y with 
our process \/N\E^ + y/N^E^ . The further identifications A = T, N xT, N , 
B = 0, / = In Z immediately entail that Hypothesis (1) of 0] reduces to 
( |TT| ) and Assertion (3) to our formula fl3"6|) because Hypothesis (2) is just 
convexity of In Z. 
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3 Examples 



3.1 The SK and even p-spin models 

For the sake of completeness we recover here the Guerra-Toninelli result 



GuTo]. First note that by the definition we have 



q N {(7,r) - -^rgjVi(7Ti(o"),7ri(T)) - ^9v 2 (7T2(o r ),vr 2 (r)) = 0. (37) 



so that ( |T7D holds as an equality for p = 1 (the random field model). By 
this means that the random field model free energy density doesn't depend 
on the size: «at = a±. For p = 2u (SK corresponds to u — 1) formula fl3~7| ) 
together with the convexity of the function x — > x 2u implies (p"7|): 

<&(',t) - ^^Ma)Mr)) ~ (*), ^(r)) < . (38) 

For the standard p-spin model defined as 



Ea = \l 2^ J h,.,i P °h ( 39 ) 



n<...<i 



p 



we refer to fCuTo 



3.2 The REM 

The model is defined by: 



kv{E a E a ,) = 5 a y. (40) 



Condition (IT7D is verified because it becomes 



^a,a' < -JT^^^ia') + -J^^ 2 (a),n 2 (s') • (41) 
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In fact if a = a' the previous formula is an identity. If a ^ a' the left hand 
side is but the right hand side is not always zero. Let us take for instance 
a = (+,+) and a' = (+,-), 7n(+,+) = +, 7Ti(+-) = +, tt 2 (+,+) = +, 
7T2(+, — ) = — . In that case the left hand side is zero and the right hand side 
is 1/2. 



3.3 The GREM 



To show the inclusion in our scheme of the Derrida- Gardner GREM [peGa 
let us first recall its construction. The GREM considers 2^ Gaussian random 
energies H(n) = ^/NE^. Their covariance is specified after the assignment of 
a rooted tree with n layers and 2 N leaves, n < N. The root furcates into 
branches, the vertices at the end of the first layer furcate into branches 
etc., up to the vertices at the end of the n — 1 layer which a^- furcate into 
the 2^ leaves. 

Remark 7 The topological constraint over the successive furcations which 
end up on 2^ leaves implies YHi=i a f = ^ ■ Each af is an integer which by 
the previous formula divides 2 N . By the fundamental theorem of arithmetics 
otf = 2 ki . Here hi, i = 1, . . . , n is a non-negative integer, and ki + k 2 + ...k n = 
N. In other words: given any tree with 2 N leaves the construction allows only 
for furcations in powers of 2 at each layer. The cofficients must depend 
on iV: in fact, cti = 2n and the only iV-independent choice of the vector a is 
obtained for ki = j- where the integers U have to divide iV for all N . Hence 
they must fulfill the constraint Y^h=\ t = ^ which is impossible. 
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The previous remark allows us to associate to each leave \i a spin configura- 
tion {<7i, (72, ctjv}. This can be done observing that the an(N) N = 2 kl{N) 
branches emerging from the root identify canonically the configurations of k\ 
spins, the successive branches the configuration of k 2 spins and so on. We 
have in this way associated to each leaf either a path (the only one joining 
the root to it) or a spin configuration. The model is finally specified by the 
formula E(fj) = J27=i where the e« are thrown according to n Gaussians 
with Av(ej) = and Av[(ej) 2 ] = af to each branch of the tree we associate 

an independent e whose distribution depends (through its variance) only at 

1-1 

which layer starts the branch. Defining = ^ a« , (v^ = and = 1) 

i=i 

it is immediate to prove that if two paths /i and v merge at the level I we 
have Ay (E^En) = v^. For fixed n and iV this construction is exactly the 
Derrida- Gardner process over a tree T„ jN ] we will denote it {£ ,T nN }. 

Theorem 1 entails existence of the thermodynamical limit for the GREM, 
in the sense that if {£, T Ui n} is assigned for a given n and all N > n, and the 
sequence of {T n ^} is increasing i.e. fcj(iV) > fcj(M) for N > M then its free 
energy density is (at fixed n) decreasing (and bounded) in N. To show this 
assertion, starting from a process {£, Tn ^} we build the process {£$ , T Hi n} 
with N = Ni + N 2 in the following way: at each vertex of the tree T ntNl 
sitting on the layer % we increase the multiplicity of the furcation by a factor 
2k i {N)-k i {N 1 ) ass ig n i n g the same value ep to all newly introduced branches. 
By construction the new process will enjoy the property 
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We apply the same construction to build {£n 2 , %i,n} and we have 

Av(^X 2 i)) ^ vil) ■ ( 43 ) 
It is now straightforward to verify that conditions (42) and fl43|) imply fllTI). 
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